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Correlated excitations in a plane circular configuration of identical atoms with parallel dipole moments are 
investigated. The collective energy eigenstates, which are formally identical to Frenkel excitons, can be com- 
puted together with their level shifts and decay rates by decomposing the atomic state space into carrier spaces 
for the irreducible representations of the symmetry group Zjv of the circle. It is shown that the index p of these 
representations can be used as a quantum number analogously to the orbital angular momentum quantum num- 
ber / in hydrogen-like systems. Just as the hydrogen .s-states are the only electronic wave functions which can 
occupy the central region of the Coulomb potential, the quasi-particle corresponding to a collective excitation of 
the atoms in the circle can occupy the central atom only for vanishing In quantum number p. If a central atom 
is present, the p = state splits into two and shows level-crossing at certain radii; in the regions between these 
radii, damped quantum beats between two "extreme" p = configurations occur. The physical mechanisms 
behind super- and subradiance at a given radius are discussed. It is shown that, beyond a certain critical number 
of atoms in the circle, the lifetime of the maximally subradiant state increases exponentially with the number of 
atoms in the configuration, making the system a natural candidate for a single-photon trap. 

PACS numbers: 42.50.Fx, 32.80.-t, 33.80.-b 



I. INTRODUCTION 

When a collection of identical atoms is located such that 
their mutual distances are comparable to the wavelength of an 
atomic transition, the mode structure of the electromagnetic 
field is altered, at least for photon states whose energy lies 
in the neighbourhood of the associated atomic transition; at 
the same time, level shifts and spontaneous emission rates of 
the states in the atomic ensemble undergo changes. Further- 
more, the atoms in the sample become entangled with each 
other as well as with the radiation field via photon exchange. 
The sample then occupies collective states whose associated 
spontaneous emission rates can be greater or smaller than the 
single-atom decay rate and are called super- or subradiant, 
accordingly. For a collection of two-level atoms, this phe- 
nomenon was discovered as early as 1954 by Dicke [1]. 

Super- and subradiance in the near-field regime has been 
studied by a number of authors. Reviews of the subject were 
given in [2, 3]. The theory of subradiance in particular was 
treated in a four-paper series in [4-7]. Wigner functions, 
squeezing properties and decoherence of collective states in 
the near-field regime have been presented in [3, 8, 9], while 
triggering of sub- and superradiant states was shown to be 
possible in [10]. Experimentally, super- and subradiance have 
been observed in [11, 12]. 

In all these investigations the so-called "small-sample ap- 
proximation" played a crucial role: Here the atoms are as- 
sumed to be so close together that they are all subject to the 
same phase of the radiation field. "Closeness" here is mea- 
sured on a length scale whose unit is the wavelength X of the 
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dominant atomic transition under consideration. The small- 
sample approximation is then manifested by using an interac- 
tion Hamiltonian which is independent of the spatial location 
of the constituents, similar to a long-wavelength approxima- 
tion. Thus, in this approximation, super- and subradiance ap- 
pears to be a near-field phenomenon. 

In this paper we perform a detailed investigation of the 
simply-excited correlated states observed in a planar circular 
configuration of atoms beyond the small-sample approxima- 
tion: That is to say, we determine the level shifts and decay 
properties of correlated states for atomic ensembles in con- 
figurations with arbitrary radius of the circle, not restricted to 
the small-sample limit. The simply-excited correlated states 
so obtained turn out to be formally identical to Frenkel ex- 
citons, e.g., in a molecular crystal. Conceptually, however, 
Frenkel excitons differ from our collective excitations in that 
the former are usually thought of as being mediated by the in- 
teratomic Coulomb interaction only [13-16], whereas, in the 
latter case, the derealization of the excitation energy occurs 
through the exchange of transverse photons and therefore pro- 
ceeds through radiative processes. Below we shall comment 
on this conceptual difference in greater detail. 

Taking advantage of the discrete rotational symmetry Zjv 
of the circle, the complex eigenvalues and eigenstates of the 
non-Hermitean channel Hamiltonian governing the dynam- 
ics in the radiationless subspace of the system can be deter- 
mined analytically by group-theoretical methods. The physi- 
cal mechanism determining super- or subradiance of a given 
eigenstate is then discussed. In particular we show that each 
of the collective eigenstates has certain ranges of the radius 
of the circle for which it can be the maximally super- or sub- 
radiant state. We demonstrate that, for an ensemble with N 
outer atoms in the circle, there are (N — 1) states which are 
insensitive to the presence of a central atom; this means that 
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the associated wavefunctions, level shifts and decay rates of 
the quasi-particle describing the collective excitation do not 
change when the atom at the center is removed, simply be- 
cause the latter is not occupied. There are only two states 
in a circular ensemble which occupy the central atom; we 
show that they are analogues of the s-state wave functions in 
hydrogen-like systems: Just as the s states carry angular mo- 
mentum quantum number I = and therefore transform un- 
der the identity representation of 5(9(3), our collective p = 
states correspond to the identity representation of the symme- 
try group ~L N \ and just as the s states are the only ones which 
are nonvanishing at the center-of-symmetry of the Coulomb 
potential in a hydrogen-like system, so are our p = states 
the only quasi-particle states which occupy the central atom at 
the center-of-symmetry of the circle. These two p = states 
for the configuration with central atom exhibit quantum beats 
between two extreme configurations, one, in which only outer 
atoms are occupied, the other, in which only the central atom 
is excited. At certain radii the two levels cross, and no quan- 
tum beats occur; in this case the population transfer between 
the two extreme configurations is completely aperiodic. For 
other radii, the beat frequency is so much smaller than the de- 
cay rate of the collective state as a whole, that the decay effec- 
tively proceeds like an aperiodic "population swap". Finally 
we show numerically that, by increasing the number of atoms 
in the circle for fixed radius, we arrive at a domain where the 
minimal spontaneous decay rate in the sample decreases ex- 
ponentially with the number of atoms in the circle. Such con- 
figurations therefore are natural candidates for single-photon 
traps. This pronounced photon trapping may even survive in 
the presence of dephasing interactions with an environment, 
or more generally, when the system is subject to a (moderate) 
loss of quantum coherence due to external noise. 



II. SPONTANEOUS DECAY AND LEVEL SHIFTS IN A 
SAMPLE OF CORRELATED ATOMS 

In this work we do not study the most general collec- 
tive states in an atomic ensemble but focus on those states 
which are coupled to one resonant photon only; we shall call 
these states simply-excited. The simply-excited uncorrelated 
states of the sample are given as \e\,g2,--- ,gN,0) = 1 1 , 0) , 
|^i,^2, ■ • • ,eN,0) = \N,0), where \g A ), \e A ) denote the 
ground and excited level of the Ath atom, and |0) is the ra- 
diative vacuum. These states are coupled to the continuum 
of one-photon states \g\,g2, • • • ,gN,ks) = \G,ks), where (ks) 
denote the wave vector and polarization of the photon. The 
uncorrelated states |A,0) are occupied by a sample of atoms 
for which the following is true: 1.) the atoms have infinite 
mutual distance; and 2.) precisely one of the N atoms is ex- 
cited. 

If we now think of adiabatically decreasing the spatial sep- 
aration of the atoms to finite distances, the atoms will begin 
to "feel" each other via their mutual Coulomb (dipole-dipole) 
interaction and their coupling to the radiation field, and the 
states |A,0) will no longer be the stationary energy eigen- 
states. Rather, we expect that superpositions of uncorrelated 



excitations will emerge, 

N 

\C)= Y. c b\B,0) , (1) 

B=l 

with coefficients eg which are obtained by diagonalization of 
a suitable effective Hamiltonian, given below. 

Inherent in our model is the assumption that the atoms tak- 
ing part in the collective excitation have mutual distances on 
the order of magnitude of an optical wavelength, i.e., several 
hundred nanometers. As a consequence, electronic wavefunc- 
tions pertaining to different atoms do not overlap, and we can 
refrain from antisymmetrization of the total electronic state 
(Slater determinant). This justifies the use of product states 
for the simply-excited uncorrelated states |A,0). Furthermore, 
it is clear that under these conditions, the migration of exci- 
tation energy between the different atoms can proceed only 
through the exchange of transverse photons, as the Coulomb 
dipole-dipole interaction will play a role only for interatomic 
distances Rab <C h. es , where X eg is the dominant wavelength 
of the Bohr transition. This feature seems to set the corre- 
lated states (1) apart from the usual notion of Frenkel excitons 
[13-16], although both are formally identical in being super- 
positions of simply-excited product states. In the case of a 
Frenkel exciton, however, it is usually assumed that the atoms 
involved in the collective excitation are located on the lattice 
sites of a crystal and therefore have mutual distances compa- 
rable to a Bohr radius. In this case, the "excitation transfer in- 
teraction" [15] responsible for the migration of excitation en- 
ergy is Coulombic, and indeed this is assumed in all standard 
accounts on Frenkel excitons (Coulomb means direct plus ex- 
change interaction, if relevant). To quote from a standard ref- 
erence ([14], p. 19), "Excitons ... are idealized elementary 
excitations in whose consideration we ignore the delay effects 
and take into account only Coulomb interaction". This points 
out an important conceptual difference between our collective 
excitations and the standard notion of a Frenkel exciton. 

We can express this difference in even greater detail: In 
the case of the "traditional" Frenkel exciton, the system un- 
der consideration is comprised of the totality of atoms in the 
crystal, interacting through the instantaneous Coulomb poten- 
tial or its multipolar approximation. On account of this in- 
teraction, the atoms become entangled with each other, but 
certainly remain unentangled with the transverse degrees of 
freedom of the radiation field, at least as long as only quasi- 
resonant atom-photon interactions (rotating-wave approxima- 
tion) are taken into account. On the other hand, the collective 
states considered in our work here are states of the total sys- 
tem atoms+radiation, and as a consequence, in any of these 
states the atoms must be regarded as being entangled not only 
with each other but also with the transverse radiative degrees 
of freedom. This entanglement, of course, is a direct conse- 
quence of the fact that the migration of excitation energy pro- 
ceeds through exchange of a transverse photon: the excited 
atom A emits a quasi-resonant photon which, in turn, gets ab- 
sorbed by atom B, thereby exciting the latter. The electrostatic 
interaction plays a role in this process only for interatomic 
distances which are small compared to an optical wavelength. 
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Due to these conceptual differences, we hesitate to call our 
collective states Frenkel excitons. 



a) 



b) 



FIG. 1: The two circular configurations of A^+l) atoms which are 
examined: In configuration a), A' outer atoms form a regular N- 
polygon with radius r, plus an atom at the center. Configuration b) is 
the same as a), but without the central element. 

We now present the details of the computation: Let N neu- 
tral atoms be labelled by A = l,2,...,N, each atom being lo- 
calized around a fixed point Ra in space, where Ra are c- 
numbers, not operators. The atoms are assumed to be identi- 
cal, infinitely heavy, and having a spatial extent on the order of 
magnitude of a Bohr radius ao- The point charges within each 
ensemble are labelled as qAa- On taking the long-wavelength 
approximation in the minimal-coupling Hamiltonian and per- 
forming the Goppert-Mayer transformation we arrive at the 
electric-dipole Hamiltonian 

H = II +1 L v ** ( 2a ) 



A a 



A <x<p 



-£d A »E ± (R A )=//o+#/ , 



(2b) 



where the unperturbed Hamiltonian Hq, given in (2a), contains 
the sum over atomic Hamiltonians Hqa including the Coulomb 
interaction V Aa p between the internal constituents of atoms 
A = 1, . . . , N, but without the Coulomb dipole-dipole interac- 
tion V%g between atoms A and B, for A ^ B; and the normally- 
ordered radiation operators. The electric-dipole interaction is 
given in (2b), where Ej_ denotes the transverse electric field 
operator. The interatomic Coulomb dipole-dipole interaction 
V^g seems to be conspicuously absent in (2); however, it is a 
feature of the Goppert-Mayer transformation to transform this 
interaction into a part of the transverse electric field, so that 
the Coulomb interaction will emerge as a factor in the level- 
shift operator, as will be seen below in formulae (8-9). As 
a consequence of being a part of the transverse electric-field 
operator, the interatomic Coulomb interaction is now fully re- 
tarded, in contrast to the instantaneous Coulomb interaction 
in the minimal-coupling Hamiltonian. 

The process we study in this work consists of the radia- 
tive decay of a simply-excited correlated atomic state, or a 
"Frenkel exciton", which is resonantly coupled to a continuum 
of one-photon states. For this reason, two-photon- or higher- 
photon-number processes can be expected to play a negligible 
role, so that it is admissible to truncate the possible quantum 
states of the radiation field to one-photon states \ks), and the 
vacuum |0). The rotating-wave approximation is now taken 
[17], Thus, amongst the admissible single -photon transitions 
\c) — * |G,fcy) and \c) — * \AB,ks), only processes of the first 



kind are taken into account. The radiationless subspace Q is 
spanned by the states |A,0), using the same symbol for the as- 
sociated projector Q = Y!a=\ |^,0) (A,0J, while P denotes the 
projector onto the subspace of one-photon states \G, ks) . 

In the present scenario we assume that the system 
atoms+radiation is closed and therefore evolves unitarily; a 
qualitative consideration of the impact of dephasing interac- 
tions will be given in the last section VII. If, at time t = 0, 
the system is in a correlated state (1), the probability of find- 
ing the radiation at t > to be still trapped in the system is 
determined by the g-space amplitudes 



(A,0\U(t,0)\c) 



N , . 

E — / 

2%i J 



dEe-* Et 



(3) 



■(A,Q\QG(E+)Q \B,0) c B 



where the g-space Green operator QG(E + )Q is given in 
terms of the non-Hermitean g-channel Hamiltonian H (E + ) 
as QG(E+)Q = 1/{E+ - #(£+)). At the initial energy E t = 
E e + (N — l)E g , where E e ,E g are the energies of the excited 
and ground state of a single atom, we have 



X(E i+ ) = QH Q Q + hA(E i )-i-r(E i ) , (4) 
where A(£,-) is the level-shift operator 



A AB (Ei) = -hrs(x) 



X — k e „RAB 



(5a) 



3 | cosX sinX cosX 



+ 



kX 2 I 



du 



e~ u (l+u + u 2 ) 



u 2 +X 2 



and r(Ej) is the decay operator 

TAB{Ei) 

Di{X) 



TD\ (k eg R A B) , 
3 / sinX cosX sinX\ 
X^ 



2 V X 



X 2 



(5b) 

(6a) 
(6b) 



in the g-space, where T = d 2 k\ s j (37ten/z) is the spontaneous 
emission rate of a single atom in a radiative vacuum and 
E eg = E e — Eg = hck eg . The level-shift function S can be ap- 
proximated by 



S(X) 



3 \ cosX sinX cosX 



+ 



nX 2 



approx — ^ 

X 2 



X 



X 



2 



1 , 1 , 1 

arctan- + l+-ln^2 



X 3 



+ 



(7) 



A plot of the functions Dy, S and 5 apP rox is given in Fig. 2. 

[Remark: In the more general case of arbitrarily aligned 
atoms there are two correlation functions Dy,Di emerging in 
the decay matrix rather than just one. This is indicated by our 
denotation. For parallely aligned atoms, the contribution of 
L>2 vanishes.] 
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A function similar to D\, but within a classical context, is 
the dyadic Green function of the electromagnetic field [18]. 

The off-diagonal elements of the level- shift matrix Aab can 
be put into the form 



1 



47ie /i R 3 



■ ARab) 



(8) 



AB 



where /(Rab) is an analytic function of Rab which tends to \ 
as Rab tends to zero. Comparison with the Coulomb dipole- 
dipole interaction between the two atoms A and B, 



y di P _ d A • d B 3 (dU • e A g) (d B • z A b) _ d 2 



(9) 



AB 



shows that the level shifts indeed contain the dipole-dipole 
Coulomb interaction between the charged ensembles A and B. 
On the other hand, the diagonal elements of Aab diverge due to 
Coulomb and radiative self-energy; to remove this deficiency 
we renormalize [19-24] the initial energy ("bare" energy) by 
redefining 




FIG. 2: Plots of the correlation functions D\ (X) and S(X) responsi- 
ble for spontaneous decay ("D") and level shifts ("S"), respectively. 
X is equal to 2n times the interatomic distance measured in units of 
the wavelength X eg of the dominant atomic transition. D\ becomes 
maximal at zero interatomic distance while 5 tends to minus infin- 
ity there, due to Coulomb- and radiative self-energy of the atomic 
dipole. Both functions tend to zero as X — > °°, expressing the fact 
that any kind of correlation must cease to exist for infinite distance. 
The inset shows a comparison between the exact shift function S(X), 
based on eq. (5b), and the approximation S a pprox, based on eq. (7). 



■ E{ + lim Maa (EA = E t 



(10) 



and subtracting the same infinity from the level-shift matrix, 
4«(£,)-> Aab(E,-)={ t MEi) ' A / B J (ID 



where Ej is now assumed to have a finite value, namely the 
value unperturbed by the presence of N — I other atoms; as a 
consequence, £, must be assumed to have been infinite in the 
first place. The only observable level-shifts are now those due 
to interatomic interactions. 

The g-channel Hamiltonian in the uncorrected basis now 
becomes 



hY 

(A,0\X(E i )\B,0)=E i t N - 



( i S + iDi ••• \ 
S + iDi i 



V 



(12a) 



(12b) 



■I 



This matrix can be diagonalized in terms of left and right 
eigenvectors [25] 



(A,0\X(Ei)\B,0) = 

nr(k eg ) 



£<A,0|C,){^-^^^}<c;|B,0> . 



(13) 



where 



(14a) 

TJm/jp , (14b) 
For the given radius r, number of atoms N, and configuration 



(a) or (b), let | C m in) be the eigenstate with the smallest decay 
rate, 



r m in < r p for all p = 1, . . . ,N 



Then the associated correlated state 

1 



\C) 



min | >-min 



I ^min) 



(15) 



(16) 



is the most stable one with respect to spontaneous decay [say- 
ing nothing about stability against environmental perturba- 
tions], and hence is a candidate for a single-photon trap. 



III. CYCLIC SYMMETRY OF THE CIRCULAR ATOMIC 
CONFIGURATIONS 

We now construct the eigenvectors \ c p ) of 9i (E{) explicitly 
by group-theoretical means, taking advantage of the fact that 
the system has a cyclic symmetry group 



{e,T,T 2 



nN-l 



} 



(17a) 
(17b) 



where is the number of outer atoms along the perimeter 
of the circle, the generator T acts on coordinate space as a 
(passive) rotation by the angle 2%/N about the symmetry axis 
of the circle, and acts on the state space unitarily by relabelling 
the outer atoms, but leaving the central atom invariant, 



T\A) = \A- 
T \z) = \z) 



1) 



(18) 



The unitary irreducible representations of Z;v are given by 
r^)=ex P (^) , A = 1,2,...,N , (19) 
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and using appropriate projection operators, the radiationless 
2-space can be decomposed into an orthogonal direct sum 



Q = ®^Z Tp of carrier spaces T p of irreducible representa- 
tions F 1 '. Since the channel Hamiltonians of both configura- 
tions are invariant under this transformation, 



(20) 



H preserves all carrier spaces T p and hence can be diago- 
nalized on each T p separately. This simplification makes it 
possible to perform the diagonalization analytically. 



IV. EIGENSPACES OF THE GENERATOR T OF THE 
SYMMETRY GROUP 

In order to perform the diagonalization we need to construct 
basis vectors of the carrier spaces T ' p by applying the standard 
projection operators [26] <P'\ 



N-l 

E ex p 

A=Q 



2ni P A \ a 
N J 



(21) 



projecting onto T p , on an arbitrary correlated state \c) = 
c z \z,0) +La=i c a|A,0), with the result 



*' \c) = 8 p0 c z \z,0) +c £ exp \A,0) 



(22) 



where |c z | 2 + \c\ 2 = 1. It follows that, for p ^ 0, the carrier 
spaces T p are one-dimensional, and are spanned by normal- 
ized basis vectors 



C. 



1 



,> s -I>p(^)m> 



(23) 



which turn out to be simple Fourier transforms of uncorre- 
cted excitations. This statement is true for both configura- 
tions (a) and (b). We emphasize again that, here, N denotes 
the number of outer atoms. In Fig. 3 we plot the real parts 
of c p A = (A,0\c p ) for some of the /^-states with N = 50 outer 
atoms. 

If p — 0, the associated eigenspace To is one-dimensional 
for configuration (b), and is spanned by 



|Co) 



44 M 



(24a) 



while, for configuration (a), To is two-dimensional, with a 
second basis vector 



\z,0) 



(24b) 



The two states in (24) are totally isotropic under the symmetry 
group Z N in the sense that they do not pick up a phase factor 
when we operate with a group element on them. Furthermore, 
they are "extreme" states, since in (24a), only outer atoms are 
occupied, with all atomic dipoles aligned parallely, as can be 
seen in the first plot in Fig. 3; while in (24b), only the central 
atom is occupied. 
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FIG. 3: The real parts $Hec^ of the components of states (24a, 23) in 
the uncorrelated basis, for the eigenstates p = 0, 1,2,23,24,25 of a 
configuration with N = 50 outer atoms. For p ^ these coefficients 
are the same for both configurations (a) and (b). 



V. DIAGONALIZATION OF THE CHANNEL 
HAMILTONIAN ON CARRIER SPACES OF THE 
SYMMETRY GROUP 



A. Configuration (b) 



We first discuss configuration (b) without central atom. In 
this case, each of the carrier spaces T p is one-dimensional, 
and is spanned by states (23, 24a). Since 9i and % pre- 
serve these carrier spaces it follows that each of the \C P ), 
p = 0,...,N — 1 is automatically a right eigenvector of 




Level degeneracy for configurations (a) and (b) 



FIG. 4: Relative level shifts A p /T and decay rates T p /T for the states 
p — 0, 1,2,3 with N = 7 outer atoms and no central atom. From 
the second plot we see that none of the p-states can generally be 
declared as super- or subradiant; rather, for each state there exist 
ranges of the radius r for which the state | Cp) is maximally super- or 
maximally subradiant. The physical reason for this dependence lies 
in the interference due to varying phase differences between spatially 
retarded radiative contributions from all but one atom, at the site of 
the given atom. 



with eigenvalue 



2A P Xp 

v P = — f^ + 'Y 



r N 

A p = -t£ S(k eg Ru) cos ( 

L A=2 



(2%p{A-\) 



I N 

r p = r\i + Y,D 1 (ke g RiA)cos 

A=2 



N 

/27tp(A-l) 



where we have used (40). It follows that 



(25a) 
(25b) 

, (25c) 
(26) 



hence some of the eigenvalues are degenerate; the exact re- 
sult is presented in Table I. Eqs. (25b) and (25c) contain the 
level shifts and decay rates, which are given as cosine trans- 
formations of the correlation functions S and D\. A plot of 
these quantities for the first four eigenvalues for N = 7 outer 
atoms is given in Fig. 4, where the approximated shift function 
Sapprox based on eq. (7) has been used. 



B. Mechanisms for super- and subradiance 

The Figures 3 and 4 give some insight into the mecha- 
nism of super- and subradiance. We first note that the system 
atoms+radiation evolves unitarily in the present treatment, and 
since the initial state was pure, the state of the total system is 
always pure. This fact entitles us to visualize the atoms in 
the sample as being radiating coherently. This coherence is 
at the heart of super- and subradiance since it means that the 
radiative contributions emitted by the individual atoms will in- 
terfere: generally speaking, a correlated state will be super- or 





N=2n+\ 


N=2n 


degenerate 


l,...,n 
2n,. . . ,«+ 1 


p= l,-,n-l 
2n-l,...,n+l 


non-degenerate 


p=0(±) 


P=0(±),« 



TABLE I: Degeneracy of the eigenvalues for quantum numbers p = 
0(±), . . . ,N — 1, for odd and even A', for configurations (a) and (b). 
Vertically stacked p values in the second row are mutually degener- 
ate. 



subradiant when the total interference between these contribu- 
tions is constructive or destructive in the far zone (measured 
in units of X eg ) of the radiation field around the configuration. 
Based on this reasoning we now give a qualitative explanation 
for the behaviour of level shifts and decay rates in Fig. 4 in 
the small-sample limit r — > 0: 

We see in Fig. 4 that the p = state is the only one which 
has a nonvanishing decay rate in this limit: All atomic dipoles 
are aligned parallely in this state [see eq. (24a) and first plot in 
Fig. 3], and vanishing interatomic distance on a length scale of 
a wavelength means that the radiation emitted coherently by 
the sample must interfere completely constructively. Hence, 
the sample radiates faster than a single atom by a factor N, 
since the decay rate is 



limr p=0 =N-T 



(27) 



as follows from eq. (25). Conversely, the suppression of spon- 
taneous decay for the states with p > is a consequence of the 
fact that the dipoles have alternating orientations, see eq. (23) 
and the second to fifth plot in Fig. 3, so that, in the small- 
sample limit r — > 0, roughly one-half of the atoms radiate in 
phase, while the other half has a phase difference of jr.; hence 



limr p>0 : 







(28) 



This behaviour is exemplified by the plots of ri,r2,r3 in 
Fig. 4. However, the destructive interference in the (p > 0) 
states is independent of how adjacent dipoles are oriented, as 
long as r< X eg ; in this limit, the decay rate should not be 
noticeably affected by rearranging the dipoles in different pat- 
terns as long as the 50:50 ratio of parallel-antiparallel dipoles 
is kept fixed. What will be affected by such a redistribution 
is the level shift, since the Coulomb interaction between the 
dipoles in the sample may change towards more attraction or 
repulsion between the atoms. By this mechanism we can ex- 
plain the divergent behaviour of the level shifts in the small- 
sample limit: The shift of the p — state always behaves like 



limA,; 



(29) 



which arises from the Coulomb repulsion of the parallely 
aligned dipoles. On the other hand, the level shifts for the 
(p > 0) states depend on the relative number of parallely 
aligned dipoles in the immediate neighbourhood of a given 
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dipole, or conversely, on the degree of balancing the Coulomb 
repulsion by optimal pairing of antiparallel dipoles. As a con- 
sequence, states for which p is close to zero always have pos- 
itive level shifts, since the Coulomb interaction between ad- 
jacent dipoles tends to be repulsive, as seen in the first three 
plots in Fig. 3, and the behaviour of Ao, Ai in Fig. 4. On the 
other hand, states for which p is close to N/2 tend to have 
antiparallel orientation between adjacent dipoles, hence the 
Coulomb interaction is now largely attractive, which explains 
why these states have negative level shifts in the limit r <C A,. 
This is seen in the last three plots in Fig. 3 and the behaviour 
of A2, A3 in Fig. 4. 

The same mechanism clearly also governs the super- or sub- 
radiance of the sample with finite interatomic distance. In 
this case the information about the orientation of surrounding 
dipoles at sites Ra is contained in the transverse electric field, 
which arrives at the site Ri with a retardation \Ra — Ri|/c. 
Thus, in addition to the phase difference imparted by the co- 
efficients c p A , there is another contribution to the phase from 
the spatial retardation, which accounts for the dependence of 
the level shifts and decay rates on the radius r. However, 
apart from this additional complication, the physical mecha- 
nism determining whether a given state is super- or subradiant 
is clearly the same as in the small-sample limit, and can be 
traced back to the mutual interference of the radiation emitted 
by each atom, arriving at a given site Ri . 



C. Configuration (a) 

Now we turn to compute the complex energy eigenvalues 
and eigenvectors for configuration (a) with central atom. For 
p = 1 , . . . , N — 1, the basis vectors carrying irreducible repre- 
sentations of the symmetry group are the same as before, and 
are given in eq. (23). Consequently, they are also right eigen- 
vectors of the matrix ; the associated eigenvalues p p turn out 
to be the same as for configuration (b) and thus are given by 
formulas (25). This result means that the presence or absence 
of the central atom is not "felt" by the system in the modes 
I C p ), p = l,...,N—l, since the central atom is not occupied 
in these states. 

On the other hand, the eigenspace To of T corresponding 
to the p — representation is now two-dimensional and is 
spanned by r° -basis vectors (24a) and (24b). We now need 
to diagonalize the channel Hamiltonian on this subspace: the 
diagonalization of the matrix in the basis | Co) , | C z ) yields 



where 



with 



c = 



|Co+) = cosG |Co)+sin9 \c z ) , 
I Co-) = -sin9 I Co) + cos 9 \c z ) , 

2 ; — c 
2VNM(k eg r) 



(30a) 
(30b) 

(31) 



lL 2 M(V? 1A ) + yj [Z N A=2 M(k eg R lA )\ 2 +4NMHk eg r) ' 

(32) 



and M — S + iD\. The eigenvectors (30) are not orthogonal, 
consistent with the fact that the matrix 3(_ is not Hermitean. 
However, (30) form an orthonormal system together with the 
left eigenvectors 



(c *+| = cos9 (Co| + sin§ (c z \ 
(C *_| = -sin§ (Co|+ cos6 (c z \ 

The associated eigenvalues po± (of ?0 are 
1 N 

PQ±=i+~Y, M {kegRlA)± 



(33) 



' A=2 



2 \IL 



A=2 



(34) 



+ 4NM 2 (k eg r) 



where r is the radius of the circle. From (34) we obtain the 
level shifts Ao± and decay rates To± of the p = states as in 
eq. (25). 

The definition of eigenvalues p.o± as given in eq. (34) is not 
yet the final one, however. In (34) we have assigned po + to 
the square root with positive real part. From (14a) we see that 
the po+ so defined always comes with a negative level shift, 
and therefore the associated real energy Ej + Mo+ is always 
smaller than the energy associated with /jo-- This state of af- 
fairs would be acceptable as long as the levels would never 
cross; but crossing they do, as can be seen in Fig. 5: For radii 





- A o + 




— A o- 




-- A o = 

















- r o + 






/ \ 


















--- r o- 
























-- r o ' 




/ N 























0.5 1 



1.5 



2 2.5 

r [X ] 
L eg J 



3.5 



FIG. 5: Plotted are relative quantities A ( ± j/r and T ( ± )/r. For 

N = 10 outer atoms, the level shifts Aq± exhibit level crossing twice 
on every wavelength, for r > X eg . At the radius r of a crossing, the 
beat frequency (37) vanishes, preventing any oscillatory population 
transfer between | Co) and |c z ). The states |Co±) are defined such 
that, as r — + 0, the level shifts Ao± tend to ±°°, respectively. — 
For comparison, level shifts and decay rates of the state | Co) without 
central atom are included. 

greater than a certain lower bound ro, which depends on the 
number N of outer atoms, we see two level-crossings per unit 
wavelength. For 7Y = 10 this radius is roughly ro ~ A,. At 
each crossing we have to reverse the assignment of square 
roots to eigenvalues in order to obtain smooth eigenvalues. 
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We therefore have to redefine /jq±, and the associated eigen- 
vectors, in order to take account of this reversal. A look at 
Fig. 5 shows that, for r < 0.7 X, no crossings occur, so that 
we can uniquely identify the eigenvalues by their behaviour 
in the limit r —* 0. Thus, we finally define /iq± to be that 
eigenvalue whose associated level shift Ao± = —TfRe/jo±/2 
tends towards ±°°, respectively. Physically, the dipoles in the 
state |Co+) /||Co+| are aligned parallely, which in the limit 
r — > produces Coulomb repulsion, and hence the positive 
level shift. It follows that the state (p =0+) is the natu- 
ral analogue of the (p = 0) state in configuration (b), since 
they both have the same behaviour at small radii. It is also 
expected to decay much faster than a single atom, an expec- 
tation which is confirmed by the behaviour of the imaginary 
part Jm/Jo+ = To+ /T, which tends to ~ N + 1 as r — ► 0. Again 
this follows the pattern of the p = state in configuration (b). 
A visual comparison between states (p = 0+) and (p = 0) is 
given in Fig. 5 for A' = 10 outer atoms. 

On the other hand, in the state | Co—) /||Co-||, the central 
atom is now oriented antiparallelly to the common orientation 
of the outer dipoles, and, in the limit r — > 0, is much stronger 
occupied than the outer atoms, as follows from eq. (30b). 
Hence, as a consequence of Coulomb attraction between the 
outer atoms and the central atom, the energy is shifted towards 
— °°, and at the same time, the system has become extremely 
stable against spontaneous decay: This is reflected in the fact 
that, as r —* 0, the decay rate To- tends to zero as well. 

As mentioned above, the states with higher quantum num- 
bers p = 1,. . . ,N — 1 are the same as in configuration (b), 
and have the same eigenvalues. The two levels p = (0±) are 
non-degenerate, except for accidental degeneracy, and also are 
non-degenerate with the (p > 0) levels. As a consequence, the 
level degeneracy is similar to case (b), and is again expressed 
in Table I. 



D. Quantum beat between unperturbed p = states 

As follows from eqs. (30), the true eigenstates \ Cq±) are in 
general linear combinations of the "unperturbed" irreducible 
basis vectors | Co), eq. (24a), and \c z ), eq. (24b); as a conse- 
quence, the true eigenstates will give rise to a quantum beat 
between |Co) and \c z ). In order to determine the beat fre- 
quency we compute the probability of finding the system at 
time t in the extreme state | Co) if initially it was in the other 
extreme state \c z ), 



The last line shows that the oscillation between the two states 
occurs at the beat frequency 



A 0+ - A - 



(36) 



The beat frequency can be expressed in terms of the function 
M = S + iD 1 as 



cor = ^9te , 



-i 2 



\ £M(£ eg /?i A ) 

N L A=2 



4NM 2 (k eg r) . (37) 



For a given number of outer atoms, there can exist discrete 
1.5 . 

N=10 
N=60 
N=100 




FIG. 6: The beat frequency eq. (37) for three different values N = 
10,60, 100 of outer atoms, for radii between 2 and 5 wavelengths. 
At the radius of a level crossing Ao+ = Ao_, the beat frequency van- 
ishes. It can be seen that the greater the number of outer atoms N, the 
greater becomes the radius ro beyond which level crossing occurs. 

radii at which the beat frequency (x)r vanishes. From for- 
mula (36) we see that this is the case precisely when the two 
levels cross, and hence the normalized states |Co+) and |Co ,-) 
are degenerate in energy. This is demonstrated in Fig. 6. 

The beat frequency (Or so computed has to be treated with 
a grain of salt, however. The reason is that the dynamics in 
the radiationless g-space does not preserve probability flux, 
since the latter decays into the P-space when occupying the 
modes of outgoing photons. This is reflected in the presence 
of damped exponentials in formula (35). Depending on the 
number of atoms involved and the radius of the circle, this 
damping may be so strong, compared to the amplitude of the 
beat oscillations, that the dynamical behaviour effectively be- 
comes aperiodic, i.e., exhibits no discernable oscillations. An 
example is given by the N = 10 plot in Figure 7. 



E. Analogy between p = states and hydrogen-like s states 



P\c z )^\co)(t) = (Co\U(t,0)\c z ) 



sin 9 cos 9 



X COS 



(A 0+ -A _)i 



(35) 



It is interesting to note that for the p ^ states, the central 
atom is unoccupied, irrespective of the radius of the circle, or 
the number of atoms in the configuration. This means that 
the central atom takes part in the dynamics only in a p = 
state. This is strongly reminiscent of the behaviour of scalar 
single-particle wavefunctions in a Coulomb potential, such as 
a spinless electron in a hydrogen atom: In this case, the elec- 
tronic wavefunction vanishes at the origin of the coordinate 
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- N=10 






- N=30 






N=50 











t [r- 1 ] 

FIG. 7: The probability P(t)/\ sin9cos9| 2 for the transition \c z ) 
| Co), as given in eq. (35), for N = 10,30,50 outer atoms, at a ra- 
dius r = 0.45 X eg . The beat oscillations for N = 10 are practically 
invisible, making the transfer effectively aperiodic. 



system, i.e., at the center-of-symmetry of the potential, for 
all states with orbital angular momentum quantum number I 
greater than zero. On the other hand, in the case of our planar 
atomic system, the circular configurations also have a center- 
of-symmetry, namely the center of the circle. We can interpret 
the "Frenkel excitons" \c p ) / \ \C P \\ as states of a single quasi- 
particle, which is distributed over the set of discrete locations 
R z ,Ra corresponding to the sites of the atoms. Then the am- 
plitudes (A,0| C p ) — c p A play a role analogous to a spatial wave- 
function (x|\|/) = \|/(x); and just as the hydrogen-like wave- 
functions vanish at the origin for angular momentum quan- 
tum numbers / ^ [27], so vanish our quasi-particle wave- 
functions at the central atom for all quantum numbers other 
than p — 0. In both cases, the associated wave functions are 
isotropic: The s-states transform under the identity represen- 
tation (/ = 0) of SO(3) in the case of hydrogen-like systems, 
and states | Co(±)) under the identity representation (p = 0) of 
Zjv in the case of our circular configurations. This means that 
the quantum number p is analogous to the angular momen- 
tum quantum number I in the central-potential problem; with 
hindsight, this may not surprising, since both quantum num- 
bers p and / are indices which label the unitary irreducible 
representations of rotational symmetry groups 1 N and SO(3). 



We also see that the slope is a function of the radius r. 

We note that the four curves in Fig. 8 imply the existence of 
a common critical interatomic distance: For large N, the next- 
neighbour distance R nn between two atoms on the perimeter 
of the circle is roughly equal to 2nr/N; if we compute this 
distance for the pairs of values (r,N) as found above we ob- 
tain R nn = 0.45,0.47,0.48,0.48 X, respectively. We see that a 
critical next-neighbour distance of R c « 0.5 X presents itself: 
If, for fixed radius r, the number of atoms in the configura- 
tion is increased, the next-neighbour distance R nn decreases; 
as soon as R n „ = R c is reached, the order of magnitude by 
which spontaneous emission from the maximally subradiant 
state is suppressed becomes approximately proportional to N. 
Based upon this reasoning we see that, for a given radius r, 
the critical number of atoms in configuration (b) is given by 
N =^Y. Then, we infer from Fig. 8 that, approximately, 

r min ~r-e--< r ^ N -^ , foiN>N , (38) 

where s(r) determines the slopes of the curves in Fig. 8; this 
function decreases monotonically with r. We must have the 
limit s(r) °°> 0, because for large radii all correlations be- 
tween atoms must cease to exist, and hence r m ; n — > F in this 
limit. Let x = 1/T denote the lifetime of the excited level in a 
single two-level atom; then formula (38) tells us that the life- 
time X m ; n of a maximally subradiant state pertaining to r m i n 
increases exponentially with the number of atoms, 

Tmin^V^-^ , forN>N . (39) 

Theoretically, at least, we can therefore keep a photon trapped 
in a circular configuration for arbitrarily large lifetimes, just 
by increasing the number of atoms in the circle, hence de- 
creasing the next-neighbour distance. In practice, of course, 
a lower limit would be attained when the wavefunctions of 
adjacent atoms start to overlap. Furthermore, the impact of 
the environment on the coherence of the maximally subradi- 
ant state must be taken into account: 



VI. PHOTON TRAPPING IN MAXIMALLY SUBRADIANT 
STATES 

In this section we are interested in the photon-trapping ca- 
pability of maximally subradiant states, for large numbers of 
atoms in the circle. To this end we choose a fixed radius, 
increase the number of atoms in the configuration gradually, 
and, for each number N, compute the decay rate r m ; n of the 
maximally subradiant state for the given pair (r,N), in con- 
figuration (b) only. We then expect a more or less monotonic 
decrease of r m i n as N increases. But what precisely is the 
law governing this decrease? A numerical investigation gives 
the following result: In Fig. 8 we plot the negative logarithm 
— ln(r m j n /r) of the minimal relative decay rate at the radii 
r = 1, 1.5,2,2.5 X for increasing numbers of atoms. We see 
that from a certain number N = N onwards, which depends 
on the radius, — ln(r m i n /T) increases approximately linearly 
with N; in the figure, we have roughly N ~ 14 for r = X, 
N = 20 for r = 1 .5 X, N = 26 for r = 2 X, N = 33 for r = 2.5 X. 



- - r=1 




■ --- r=1.5 




r=2 




- r=2.5 







10 20 30 40 

N 



FIG. 8: The negative logarithm of the minimal relative decay rate, 
being proportional to the order of magnitude of suppression of spon- 
taneous decay, as a function of the number of atoms N, for fixed radii 
as given in the legend. We see that for each radius there exists a criti- 
cal number N of atoms in the configuration beyond which the order of 
magnitude by which spontaneous emission is suppressed is roughly 
proportional to N. In the Figure, these critical values are roughly at 
N = 14,20,26,33 atoms. Beyond these values, the next-neighbour 
distance, which is roughly equal to 2nr/N, becomes smaller than the 
critical distance X eg /2. 
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VII. DEPHASING: THE IMPACT OF AN ENVIRONMENT 

The results presented so far have been idealized in that it 
was assumed that the system was closed, being protected from 
random noise and loss of quantum coherence into the envi- 
ronment and thus evolving unitarily. In practice, of course, 
this cannot be the case any longer. Through interactions and 
subsequent entanglement with an environment, the subsystem 
atoms+radiation will loose the purity of the initial state, which 
will turn into a mixture. This loss is associated with the intru- 
sion of random elements into the evolution of the subsystem, 
of which two generic types may be conceptualized: 

Firstly, the sites Ra,R z of the atomic dipoles will undergo 
random fluctuations at finite temperature; these fluctuations 
will impart small random phases on the radiative contributions 
emitted from the individual atoms in the configuration, so that 
the radiation from the sample will no longer be completely 
coherent. This reduction in coherence will reduce both the 
super- and the subradiance of collective states, since it is the 
interference based on coherence which is responsible for these 
effects. We therefore expect that the radiative decay rates of 
the superradiant states become smaller, while those of the sub- 
radiant states will increase. Nevertheless, this effect may be 
expected to be small, since the amplitudes of oscillations of 
the atomic sites will be on the order of a Bohr radius, while, 
in the present scenario, the wavelengths of the electronic tran- 
sitions are in the optical regime. 

Secondly, quantum coherence of the subsystem 
atoms+radiation will be reduced by phonon-induced or 
collisionally-induced dephasing in the phases of the quantum 
states, produced by small random shifts in the unperturbed 
ground- and excited-state energies of the single atoms due 
to the coupling to the environment. It is possible to take 
these random elements into account phenomenologically by 
averaging over the acquired random phase shifts; in the limit 
of infinitely short memory (Markoff approximation) we then 
arrive at a dephasing-decay rate which must be added to the 
radiative decay rate [28]. 

In such a case the evolution of the subsystem 
atoms+radiation must be described by an appropriate 
master equation for the associated reduced density operator. 
For example, in a solution, the typical nonradiative decay 
channel will be comprised of conversions of the exciton 
energy into thermal vibrational energy, in other words, a 
transition from the electronic excitation into a large number 
of phonons. The phase of the coherent, wavelike motion of 
the exciton, as exemplified in Fig. 3, then decays on account 
of the interaction with intra- and intermolecular vibrations. 
These vibrations, in turn, give rise to fluctuations of the 
excitation energies and of the interaction matrix elements 
and must be taken into account in the Hamiltonian by 
an additional, stochastically time-dependent contribution. 
Averaging the equation of motion of the density operator over 
the fluctuations results in the stochastic Liouville equation 
describing the exciton dynamics [29-35]. The phonons can 
also be treated as a quantum-mechanical heat bath [36-39]. 

In each case additional, nonradiative decay channels open 
up, with associated decay rates. Only when these decay rates 



are small compared to the idealized radiative decay rates stud- 
ied in this work will the results obtained here continue to make 
sense. We hope to be able to come back to this point in the 
future. 



VIII. SUMMARY 

Collective excitations together with associated level shifts 
and decay rates in a planar circular configuration of identical 
two-level atoms with parallel atomic dipoles are examined. 
The relation between these states to traditional Frenkel exci- 
tons is discussed. The state space of the atomic system can 
be decomposed into carrier spaces pertaining to the various 
irreducible representations of the symmetry group Zw of the 
system. Accordingly, the channel Hamiltonian on the radia- 
tionless subspace of the system can be diagonalized on each 
carrier space separately, making an analytic computation of 
eigenvectors and eigenvalues feasable. Each eigenvector can 
be uniquely labelled by the index p of the associated represen- 
tation. For quantum numbers p > the circular configuration 
is insensitive to the presence or absence of a central atom, so 
that the wavefunction of the associated quasi-particle describ- 
ing the collective excitation occupies the central atom only in 
a p = state. It is explained how this feature is analogous 
to the behaviour of hydrogen-like s-states in a Coulomb po- 
tential. The presence of a central atom causes level splitting 
and -crossing of the p = state, in which case damped os- 
cillations between two "extreme" p = configurations occur. 
For strong damping, the population transfer between the two 
extreme configurations is effectively aperiodic. Finally, the 
behaviour of the minimal decay rate in maximally subradi- 
ant states for varying numbers of atoms in the configuration 
is investigated; a critical number of atoms, corresponding to a 
next-neighbour distance of X eg /2 on the circle, exists, beyond 
which the lifetime of the maximally subradiant state increases 
exponentially with the number of atoms in the configuration. 

IX. APPENDIX 

Here we prove a technical result which is used in the main 
part of the paper: 

Let Ria denote the distance between atoms 1 and A (i.e., 
outer atoms only). Let F be any function of this distance, 
F = F(R\a). Let p be an integer. Then 

£F(/? 1A )sin^(A-l))=0 . (40) 

Proof: 

The sum S can be written as 
S= Y,F(R hN -A+2)sin(^f-(N-A + lfj . (41) 
The sines are equal to 

sin(^(Ar-A + l))=-sin(^(A-l)) , (42) 
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while the distances satisfy the equations Acknowledgments 

Ri,N-a+2=Ria ■ (43) 
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